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Abstract. We study the time evolution of a Bose-Einstein condensate in an 
accelerated optical lattice. When the condensate has a narrow quasimomentum 
distribution and the optical lattice is shallow, the survival probability in the ground 
band exhibits a steplike structure. In this regime we establish a connection between 
the wave-function renormalization parameter Z and the phenomenon of resonantly 
enhanced tunneling. 

PACS numbers: 03.75.Lm, 03.65.Xp, ll.10.Gh 
1. Introduction 

Renormalization effects are among the most subtle and interesting in quantum field 
theory pQ and go to the heart of our understanding of fundamental interactions. Free 
field operators create and annihilate single particles with unit probability. The presence 
of interactions makes things much more complicated and requires powerful analytical 
tools, in order to compute the propagator and the S'-matrix through the self-energy 
function and guarantee probability conservation. In the Kallen-Lehmann representation 
[2] , probability conservation is enforced via wave function renormalization, that consists 
in a rescaling of quantum fields by a factor \/~Z to take into account the effects of 
interactions. 

The propagator yields exponential decay (Fermi "golden" rule) due to the 
contribution of a pole on the second Riemann sheet in the complex energy plane [3] , but 
gets also a contribution of order (coupling constant) 2 from a contour integration in the 
complex energy plane, that modifies the exponential decay law both at short and long 
times, yielding the characteristic quadratic and power-law behaviors. The quantity \[Z 
is the modulus of the residue of a pole of the propagator in the second Riemann sheet 
in the complex energy plane and is different from unity when the latter is computed 



Unearthing wave-function renormalization effects in the time evolution of a Bose-Einstein condensate^ 

survival 
probability 




time 



Figure 1. Survival probability of an unstable system. The exponential decay law 
Pit) ~ Pz(t) — Ze~ 7t is only valid at intermediate times and its extrapolation at 
t = yields wave function renormalization Z. In general Z =/= 1. 

beyond leading order in the coupling constant. For a stable state one has Z < 1 (due 
to probability conservation in the Kallen-Lehmann representation), but for an unstable 
state Z is unconstrained and can be > 1 jl]. This is a common feature of quantum decay 
processes arising from the coupling of a discrete level to a continuum. In such a case 
the discrete level moves away from the real axis and acquires an imaginary part, which 
accounts for its decay rate. Since Z is the square of the overlap between the initial state 
and this non-normalizable complex-energy vector, known as "Gamow state" [3 El [7] , 
one can have Z > 1. 

A typical time evolution of the survival probability of an unstable system in its 
initial state is displayed in Fig. [TJ The exponential law is approximately valid in 
an intermediate time region. It is preceded by a quadratic decay (Zeno region) and 
superseded by a power law. Notice that the extrapolation of the exponential to t — 
yields a value (wave function renormalization Z) that is in general 7^ 1. We shall write 

P(t) = Ze~ Jt + additional contributions, (1) 

where the additional contributions are typically second order in the coupling (to a 
continuum of states to which the system decays) and are important at short and long 
times. At intermediate times the exponential is assumed to dominate and characterize 
the evolution 

Pit) ~ P z {t) = Ze-*, (2) 

with 7 given by the Fermi golden rule. 

The short-time quadratic (and hence non-exponential) decay was first experimen- 
tally observed by the group led by Mark Raizen in Texas [8] [9] . In this pioneering work 
focus was on the occurrence of the quantum Zeno effect and its inverse [TQI E] due to 
repeated measurements. In this article we shall focus on the same physical system used 
in [8j [9], but under very different physical conditions, both in terms of initial state and 
physical parameters [12]. This will enable us to unearth wave-function renormalization 
effects. 
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Figure 2. Survival probability P(t) of the wave packet of the atomic cloud in the 
lowest energy band of the accelerated optical lattice. The numerical solution of the 
Schrodinger equation is the full line. The experimental data (open circles) are the same 
as in Ref. [T2j and the (red) dashed line is the exponential fit based on Pz(t) in Eq. 
(35 ), whose crossing with the y axis yields the value of Z. The slope of the exponential 

In this plot Vq 



decay yields the decay rate 7 (Fermi Golden rule) 
yielding Z < 1. 



1-EreC) Fq — 0.383, 



Let us briefly sketch the main ideas to be explored in this note. We shall investigate 
the time evolution of a wave packet (a "cloud" of atoms in a Bose-Einstein condensate) 
in a tilted optical lattice. The wave packet, initially in the lowest energy band, is narrow 
in momentum space and feels a constant acceleration due to the lattice tilt. When it 
reaches the edge of the band, it can make a Landau-Zener transition [13] to the upper 
band. The survival probability in the lowest energy band is therefore characterized by 
plateaus (when the packet evolves in the band) and sharp transition regions (when the 
packet loses probability to the upper band). A typical example is displayed in Fig. [5J 
One clearly sees that in practice, the dynamics can be rather different from the idealized 
situation depicted in Fig. [Tj The overall process is involved and the quantum dynamics 
depends on a plethora of factors, ranging from probability losses to phase coherence 
effects. We shall see that Landau-Zener transitions [13] and Stiickelberg's oscillations 
[H] are two facets of this complex evolution and that by scrutinizing the features of the 
survival probability of the wave function that collectively describes a cold atomic cloud, 
one can extract crucial information on wave-function renormalization effects. 

In this article we shall endeavor to give an elementary introduction to this problem 
and discuss in more detail the model adopted in Ref. [12] to study these effects. We 
shall also scrutinize the convergence of the decay rate and wave-function renormalization 
parameter to their asymptotic values. The main idea will be to use Landau-Zener 
transitions as a benchmark for the study of wave-function renormalization effects. 
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2. Landau-Zener tunneling 



A Landau-Zener (LZ) transition takes place in a two-level system with a time-dependent 
Hamiltonian H^z{t) whose spectrum, as a function of time, is characterized by the 
presence of an avoided crossing [131 HU [T^l US] ■ In the time- independent, diabatic basis 
{|1), 1 2)}, the expectation values of the Hamiltonian are linear in time 

(l|# LZ (i)|l) = -(2|#Lz(f)|2) = -oct, (a > 0) (3) 

while the off-diagonal matrix elements are constant in time and can be assumed to be 
real by properly adjusting the phase of the basis states 

<l|#Lz(t)|2> = (2|#lz(*)|1> = (4) 
The state \ip{t)) of the system is expanded as 

|^(t)) = a 1 (t)|l)+a 2 (t)|2) (5) 
and the Schrodinger equation reads 

•4 ( a :) = (-f 1 1 1 i . (6) 




dt \ a 2 I \ 5 

The time-dependent eigenbasis of H LZ is called adiabatic basis and can be obtained by 
a rotation 



(7) 



|l(t)> \ _ / cos9(t) -sin^(t) \ / |1) 
|2(*)> J \ sin0(t) cos6(t) J y |2) 

with 6(t) = arctan(at/5)/2. The adiabatic vectors satisfy 

H LZ (t)\i(t)) = -n(t)|i(t)>, H hZ (t)\2(t)) = m\m), (8) 

with 

Q(t) = ^(atf + 5 2 . (9) 

The degeneration of the mean energy of the diabatic states at t = is reflected in an 
avoided energy level crossing in the adiabatic basis, with the two states reaching the 
minimum energy distance at t = 0. By setting 

mt)) = b 1 (t)\l(t)} + b 2 (t)\2(t)} } (10) 



from 



we get 



= -emit)), j t W)) = +m\i(t)) (ii) 




-Q{t) -ih6{t) ] h 

m{t) n(t) j [ h 



;i2) 



Therefore, in the adiabatic basis, the evolution of the system is governed by the 
Hamiltonian 

H«(t) = -n(t)a z + T(t)a y , T(t) = W{t) = (13) 
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with <j y and a z the Pauli matrices. It is clear that the nature of the coupling in the two 
bases is very different: the coupling between diabatic states is constant in time, while 
the adiabatic states are significantly coupled only near t = 0, the off-diagonal terms in 
(13) being inversely proportional to the square of the distance between the adiabatic 
levels 2Q(t). This enables one to simplify the description of the LZ transition in the 
adiabatic basis. As a useful approximation, we can replace the Hamiltonian ( 13 ) in an 
interval of width Tlz ~ 5 /a around t = with one with constant coefficients 

ha 
25' 



tfad(t) ^ #ad(0) = -6<T K + -^Vy, (\t\ < T LZ /2) (14) 



and assume that outside this time interval the adiabatic states evolve uncoupled, 
Had(t) ~ -ata z . 

The afore-mentioned transition time Tlz can be fixed by imposing that the 
probability that the system, prepared in state \l(t)) at t — > -co, evolves into |2(t)) 
at t — > +00, be given by the Landau-Zener transition probability 

P LZ = exp (-_ . (15) 



Moreover, the o z term in the Hamiltonian (14) can be neglected with good 
approximation provided 

be sufficiently small. Within this approximation, the (unitary) evolution in the interval 
(-T L z/2,T LZ /2) is governed by 



TT — \ - cos(haT LZ /25) - sin (haT LZ /25) i 

-ex Pl i £ 2 y J ~ { sm(haT LZ /25) cos (tk*T LZ /25) 1 1 ' 



with 



sin 2 ( =T LZ I = Plz- (1* 



ha r 
25 

In the next section we will discuss how to apply this model for LZ transitions to the 
problem of interband transitions in a sinusoidal lattice. 



3. Interband tunneling in a lattice 

We are interested in describing the tunneling process from the first to the second band 
of a Bose-Einstein condensate trapped in an optical lattice. We assume that almost 
all atoms are in the condensate, so that the system is described by a single-particle 
wave function ip(x,t) [TT| . Moreover, let the condensate be dilute enough so that the 
interaction between particles can be neglected. This implies that the wave function of 
the system obeys a linear Schrodinger equation. 

The experimental condition is that of an accelerating one-dimensional optical 
lattice, with constant acceleration a [18J. In the rest frame of the lattice, a particle 
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of mass m in the lattice is subjected to an external force F = ma and a potential 
iy/2) cos(27rx/(i L ), V being the lattice depth and d L the lattice period (half wavelength 
of the counterpropagating laser beams). The dynamics of the system depends on two 
dimensionless parameters [19], related to lattice depth and external force: 

F ^ with E -=U£)- (19) 

In the adiabatic approximation, no transition takes place between bands (single-band 
approximation). This is consistent if Fd L < V, namely if F < V in Eq. ( Jl9| , and leads 
us to the two-level approximation outlined in the previous section. 

We are interested in experimental setups in which the initial state is highly 
peaked around a single quasimomentum value k , that is, the width of the initial 
quasimomentum distribution is much smaller than the width 27r/di of the first Brillouin 
zone B. Under this condition, in the adiabatic single-band approximation, the average 
quasi-momentum evolves semiclassically [TBI 120] , so that at time t 

Ft 

k(t) = k + —, (20) 
n 

with negligible spread in the distribution occurring during the evolution. This yields 
Bloch oscillations in a tilted lattice with a Bloch period 

T _ 2lTh _ h 2n (21) 
Fd L E Tec F 

The initial state analyzed here has a well defined initial momentum (in £>), but can be 
distributed among different bands. At the end of each Bloch period, the amplitude in 
band a acquires the following phase with respect to the amplitude in band (3 

2tt 

A<p a/l = -—(E a (k)-Ef } (k)), (22) 

where (. . .) denotes the average over B and E^(k) is the energy of the state with 
quasimomentum k in band 7 in units E rec . The relative phases acquired by states in 
different bands are at the origin of Stiickelberg oscillations in the interband transition 
rates [T3] , 

If the initial quasimomentum distribution is very peaked around its central value, 
the interband transition can be analyzed as a Landau-Zener tunneling, since in suitable 
parameter ranges the transition from the first to the second band is concentrated at the 
edges of the Brillouin zone, where an avoided crossing occurs. In the Pisa experiment 
analyzed here p2] , this leads to an alternation of plateaus and steep transition regions of 
the survival probability in the first band (see Fig. [2]). The diabatic basis is represented by 
the momentum eigenstates, crossing at hrr/di and coupled with strength V/A, while the 
adiabatic basis is represented by the quasimomentum eigenstates. The LZ parameters 
in Eq. (15) are 

71 hF r V 

tt = w? ' = 7- (23) 
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yielding the transition probability (18) 



^o^o)=exp(-^|. (24) 



Even if the essential features of the transition are included in Eq. (24), discrepancies can 
arise between the idealized case and the real physical situation. Indeed, the periodicity 
of the lattice implies that the afore-mentioned process occurs in a finite time, and that 
in the initial and final states the adiabatic levels are not infinitely separated. The 
corrections to the LZ transition probability due to the finite duration of the process are 



discussed in [201 [19]. Other corrections to Eq. (24) should be considered if the lattice is 
not shallow. In this case, couplings to higher momentum states play an important role 
and a two-level description is no longer a good approximation. 

Moreover, there is another kind of deviation from LZ, which is the main object of 



our analysis. Since Eq. (24) is obtained under the hypothesis that only one of the two 
adiabatic states is initially populated, it is no longer valid if both states are populated. 
These deviations can be relevant even if one of the initial populations is very close to 
zero, since they scale as the square root of the smallest population, as will be discussed in 
the following. These interference effects lead to resonantly enhanced tunneling (RET): 
the transition probability is enhanced by a large factor with respect to the LZ prediction 
if the energy difference Fdj^Ai between two potential wells (d^ being the lattice spacing 
and g^Az the distance between the wells) matches the average band gap of the non-tilted 
system. 

The dynamics of interband tunneling can be schematized by separately describing 
the transition at the avoided crossing and the evolution far from the edges of the Brillouin 
zone. We will assume that the transition between the first and the second band occurs 



in a time that is negligible with respect to the Bloch time defined in (21), namely 
Tlz *C T b . We are thus assuming that the evolution inside the first band is adiabatic 
for all k, except for k ~ 7t/g?£, when a transition towards the state with the same 
quasi-momentum in the second band is possible. 

This transition can be effectively described by an evolution operator corresponding 



to the one defined in Eq. (17): 

U= ( 312 ~ Pu ) . (20) 

\ Pl2 Si 2 J 

with = \/Plz(Vq, Fq) and p X2 = yl — s\ 2 . The operator U acts on the two- 
dimensional space spanned by {|1), |2)}, the states with the same quasimomentum in 



the first and the second band, respectively, that evolve according to (20) with k = 0. 



The transition from the second to the third band can be schematized as the loss of a 
fraction 1 — s| 3 in the population of the second band towards a continuum, occurring at 
the crossing around k = 0. This assumption is justified for small values of Vq, such that 
a particle in the third (or higher) band can be considered free. The survival amplitude 
s 2 3 can be determined by imposing that its square be equal to 1 — (Vq, F ), where 
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' (Vo, Fq) the LZ transition probability at the avoided crossing from the second to 
the third band 

P^o,F )=exp(-^|). (26) 

During each Bloch cycle separating two successive transitions, the relative phase between 
the second and the first band amplitudes increases by the quantity (22), that reads 

<l>(V ,F ) = ~(AE(y )), (27) 

-TO 

where (AT) is the energy difference (in units T rec ) between the second and the first 
band, averaged over B. The effects of the dynamics in a time Tb from one transition 
to the next one can thus be modelled in the basis {|1), |2)} by an effective non-unitary 
operator 



W = I ,, . (28) 




The lack of unitarity is due to the two-level approximation adopted in our analysis and 
signals the flow of population out of the relevant two-dimensional subspace. Observe 
that the global evolution of the condensate is always coherent and does not involve any 
atomic losses. 

By making use of this simplified model, we describe the time evolution in the 
following way. At t = the condensate is in the first band, with quasi-momentum close 
to k = 0. As the lattice is accelerated, the quasi-momentum increases until it reaches 
ti/oIl at t = Tb/2, where the operator U comes into play and transfers part of the 
population to the second band. The evolution from Tb/2 to 3Tb/2 is summarized by 
the application of W. Then, the second transition occurs, and part of the population 
in the second band (minus losses towards the third band) can tunnel back to the first 
band due to the action of U, giving rise to interference effects. The same steps occur in 
the subsequent transitions. 

On a time span Tb, the dynamics of the system is determined by the action of the 
non-unitary operator 

j<t> 



U = UW=\ 312 ~ Pl2S23 t ) (29) 



P12 s 12 s 23 e 

in the basis {|1), |2)}. The order of the two operations is not relevant, since W acts 
trivially on the initial state |1) before the first transition. 



4. Transient and asymptotic dynamics 

We now look at the time evolution evolution implied by the model outlined in the 
preceding section. The state of the system before the first transition is |1). Immediately 
after the n-th transition, occurring at time t = Tb(u + 1/2), the state of the system is 

\<5> n ) = U n \l). (30) 
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Let tpi and ijj 2 be the normalized non- orthogonal eigenvectors of the matrix U belonging 
to the eigenvalues e\ and e 2 , respectively. The initial state 

|l) = Cl |^ 1 ) + c 2 |^ 2 ), (31) 

will evolve at time T^(n + 1/2) into 

\^>n) = c 1 e n l \^ 1 ) + c 2 e^ 2 ). (32) 

Due to the non-unitarity of W, the two eigenvalues are smaller than unity. Let us assume 
that the eigenvalue e\ is larger in modulus than e 2 , so that it eventually dominates in 
the iteration (32). (This condition is necessarily verified for < Si 2 , s 23 < 1.) For n 
sufficiently large, the evolution reaches an asymptotic regime, in which the state after the 
n-th transition is determined only by the state after the previous one, with a transition 
rate depending on the largest eigenvalue. Since the survival probability in the first band 
can be defined as P n = |(l|$ n )| 2 , in the asymptotic regime one gets 

P n ^\e x \ 2 P n ^. (33) 

We define the asymptotic transition rate 

7 = -ln(| ei | 2 ), (34) 

and introduce an exponential function P%{t) that coincides with the survival probability 
at the center of the plateaus, at times t = nT^'- 

Fk(t) = Zexp{-jt). (35) 

This is the dashed line plotted in Fig. [2| (Note the logarithmic scale on the ordinates.) 



Observe that Eqs. (33)-(35) are valid in the asymptotic regime. Before reaching 
it, the ratio between P n and P n -\ in Eq. (28) depends on n through a time-dependent 
decay rate 7„, as in 



P 



n+l 



1 P 
1 n- 



(36) 



It is easy to prove that the succession 7 n converges to the value 7 of Eq. (34), unless 
the eigenvalues of U have equal moduli. Indeed 

2 



In 



P 



n+l 



Pn 



In 



c l e n l +1 W l ) + c 2 e n 2 +l {\\^ 2 ) 



Cl e?<l|Vi> + c 2 e£<l|V> 2 } 



(37) 



which, by definition of 7 in Eq. (34), reads 



7 n = 7 - 2 In 



1 + 



Cl(l|</>1> 



n+l 



1 + 



(3* 



ei(lh/>i> 

Thus, if |e 2 | < |ei|, the argument in the logarithm approaches one and 7„ — > 7. 
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5. Wave function renormalization 



The parameter Z in Eq. ( 35 ) is in general different from unity, due to the transient regime 
at the beginning of the evolution. It represents the extrapolation of the asymptotic 
exponential probability back at t — 0. We derive below an expression for Z. In 



the asymptotic regime, the system evolution described by Eq. (32) corresponds to an 
evolution operator applied to an initial unnormalized vector |\[>o) — c i I V^i) : 

|$ n ) ~ Cl e?|Vi> = U n ( Cl |Vi» = tn^o>- (39) 

The Z parameter, representing the extrapolation of the asymptotic behavior back to 
t — 0, can be defined as the square modulus of the projection of the fictitious initial 
vector \^o), onto the actual initial state |1) 

Z^|(l|M/o)| 2 = | Cl | 2 Kl|^)| 2 , (40) 

which corresponds to an extrapolated "survival probability" in the subspace spanned 
by |1), evaluated at the initial time. In order to gain a qualitative understanding of 
the dependence of Z on the phase difference acquired during a Bloch cycle, let us 
compare the first and second transitions. Let Pq — 1 be the initial value of the survival 



probability in the first band. After the first transition, application of Eq. (29) yields the 
survival probability 

P 1 = s 2 l2 P = e-^Po. (41) 

At the second transition, the discrepancy with the LZ prediction becomes manifest. In 
the parameter regime of small Vq we are considering, the ratio S23/S12 is very small [12] 
and we can apply a first-order approximation, yielding 

P 2 ~ (s\ 2 - 2s 2S a$ 3 cas<i>)P 1 = e-»Pi. (42) 

This clarifies that the plateaus in Fig. [2] are not equal. If the phase is <fi = 2irj, 
with integer j, the second transition is enhanced with respect to the first one. Thus, 
comparing with Eq. (27), local maxima in the transition rate are expected when 
P o (0 = 2vrj) = (AE(V ))/j. 

This is the mechanism at the origin of wave function renormalization. A backwards 
extrapolation of the second step gives a rough estimate of the Z parameter, which we 
call Z\. 

Z\ = exp(7i - 70) ~ 1 + 2s 23 f — J cos0. (43) 



Even if Eq. (43) represents a rather crude approximation, it is very useful in an 
experimental context, where only the first few steps in the Bloch cycles are accessible. 
If the survival amplitude can be measured up to the iV-th transition, the Z parameter 
can be approximated by 

JV-l 

Z N = exp (iV 7jv - 7») -> Z , ( 44 ) 

n=0 
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Figure 3. Scaling plot of Z — 1 vs. <p of Eq. (27), derived from RET experimental 
results obtained in Ref. |12j (open circles) using a narrow atomic quasi-momentum 
distribution. The full lines are the theoretical predictions for Vq = 1,2,3,4 (a smaller 
value of Vq yielding wider oscillations). 



where the convergence to a fixed asymptotic value Z is due to (38). The convergence 
to Z is typically very fast, and the first few cycles are already sufficient to obtain an 
excellent approximation. Z is displayed in Fig. [3] as a function of <j). 



6. Conclusions 



In this article we studied Landau-Zener transitions and used them as a benchtest 
for the study of wave-function renormalization effects in quantum decay processes. 
By scrutinizing the features of the survival probability of the wave function that 
collectively describes an ultra-cold atomic cloud, we consistently defined Z and extracted 
information on its behavior as a function of the experimental parameters. 

The value of Z reflects the overlap between a discrete state and a continuum of 
states (to which the discrete state decays). Pictorially, one might say that Z detects the 
overlap Z = | (^gI^/v) | 2 between the (generalized) decaying eigenfunction with complex 
energy eigenvalue, the Gamow state ^g, and the initial (physical) state ipi Q . The case 
Z < 1 does not surprise, being in accord with the familiar Kallen-Lehmann paradigm [2] 
for stable asymptotic states. The situation Z > 1 is more curious, and is a consequence 
of the fact that ipc does not live in the Hilbert space, its norm being infinite (5J El [7], 
and its overlap Z with the initial state can exceed 1. 

It is remarkable that a quantity like Z can be directly measured and that its 
deviations from unity yields directly measurable consequences. In addition, as the 
experimental parameters are varied, Z can take values that can be smaller or larger 
than unity. 
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